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Abstract. Nematic liquid crystals at rough and fluctuating interfaces are analyzed within the Frank elastic 
theory and the Landau-de Gennes theory. We study specifically interfaces that locally favor planar an- 
choring. In the first part we reconsider the phenomenon of Berreman anchoring on fixed rough surfaces, 
and derive new simple expressions for the corresponding azimuthal anchoring energy. Surprisingly, we find 
that for strongly aligning surfaces, it depends only on the geometrical surface anisotropy and the bulk 
elastic constants, and not on the precise values of the chemical surface parameters. In the second part, 
we calculate the capillary waves at nematic-isotropic interfaces. If one neglects elastic interactions, the 
capillary wave spectrum is characterized by an anisotropic interfacial tension. With elastic interactions, 
the interfacial tension, i.e., the coefficient of the leading q'^ term of the capillary wave spectrum, becomes 
isotropic. However, the elastic interactions introduce a strongly anisotropic cubic term. The amplitudes 
of capillary waves are largest in the direction perpendicular to the director. These results are in agreement 
with previous molecular dynamics simulations. 

PACS. 61.30 Hn Liquid crystals: Surface phenomena - 61.30 Dk Continuum theories of liquid crystal 
structure - 68.05 -n Liquid-liquid interfaces 



1 Introduction 

Nematic liquid crystals are fluids with long range orien- 
tational order ,1,. Compared to interfaces and surfaces in 
simple fluids, surfaces of nematic fluids have several pecu- 
liarities. First, the interface orients the nematic fluid j2| 
EE]. This phenomenon, called surface anchoring, is quite 
remarkable, because it implies that the surface has di- 
rect influence on a bulk property of the adjacent fluid. 
It also has well-known practical applications in the LCD 
(liquid crystal display) technology 5 . Surface anchoring 
is driven by energetic and geometric factors, and depends 
on the structure of the surface. Second, the oriented ne- 
matic fluid breaks the planar symmetry of the interface. 
This should influence the properties of free interfaces, e.g., 
the spectrum of capillary wave fluctuations. Third, the ne- 
matic fluid is elastic in a generalized sense, i.e., fluctua- 
tions of the local orientation (the director) have long range 
correlations [Hj . Since interfacial undulations and director 
fluctuations are coupled by means of the surface anchor- 
ing, this should introduce long range elastic interactions 
between the undulations. Hence interesting effects can be 
expected from the interplay of surface undulations and 
director fluctuations in liquid crystals 7_. 

In the present paper, we examine these phenomena in 
the framework of two continuum theories - the Frank elas- 
tic theory and the Landau-de Gennes theory. To separate 
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the different aspects of the problem, we flrst consider a ne- 
matic liquid crystal in contact with a fixed, rough or pat- 
terned surface (Sec.|21). In LCDs, alignment layers are of- 
ten prepared by coating them with polymers (polyimides) 
and then softly brushing them in the desired direction of 
alignment ■ Assuming that brushing creates grooves 

in the surface (TUI, the success of the procedure indicates 
that liquid crystals tend to align in the direction where the 
surface modulations are smallest. Similarly, liquid crystals 
exposed to surfaces with stripelike gratings were found 
to align parallel to the stripes |11II12II13II14| . Molecular 
factors of course contribute to this phenomenon, but the 
effect can already be explained on the level of the elas- 
tic theory. This was first shown by Berreman ^3], and 
the theory was later refined by different authors |16II17I 
118) . Here we reconsider the phenomenon and derive simple 
new expressions for the anchoring angle and the anchoring 
strength. 

In the second part (SecOJ, we consider the capillary 
wave spectrum of free nematic/isotropic interfaces. Cap- 
illary waves are soft mode fluctuations of fluid-fluid in- 
terfaces, that are present even in situations where fluctu- 
ations can otherwise bee neglected. They were flrst pre- 
dicted by Smoluchowski and the theory was later 
worked out by various authors |2UII21II22II23II24I|25II26| . Since 
then they were observed in various systems experimen- 
tally (27u28ll29ll!'i()ll!'{l| as well as in computer simulations p2l 

lary wave spectrum is governed by the free energy cost 
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of the interfacial area that is added due to the undu- 
lations. Assuming that the fluctuating interface position 
can be parametrized by a single- valued function h(x,y), 
and that local distortions dh/dx and dh/dy are small, 
the thermally averaged squared amplitude of fluctuations 
with wavevector q is predicted to be 

(IMq)P> = (1) 

where a is the interfacial tension. Note that (|ft,(q)p) di- 
verges in the limit g — > 0, hence the capillary waves with 
long wavelengths are predicted to be quite large. In real 
systems, however, the two coexisting fluids usually have 
different mass densities, and the gravitation introduces a 
low- wavelength cutoff in Eq. . 

In the last years, capillary waves have attracted re- 
newed interest in the context of soft condensed matter sci- 
ence. This is mainly due to the fact that typical interfacial 
tensions in soft materials are low, typical length scales are 
large, and coexisting phases often have very similar mass 
densities. Therefore, the capillary wave amplitudes in soft 
materials tend to be much larger than in simple fluids. 
For example, capillary waves were shown to have a signif- 
icant effect on experimentally measured interfacial widths 
in polymer blends |41II42II43| . Recently, Aarts et al. have 
even succeeded in visualizing capillary waves directly in a 
colloid-polymer mixture 

Liquid crystals are a particularly interesting class of 
soft materials, because of the additional aspect of orien- 
tational order. The present study is partly motivated by 
a recent simulation studies of the nematic/isotropic inter- 
face in a system of ellipsoids |2Zj, where it was found that 
(i) the capillary wave spectrum is anisotropic, and (ii) the 
interface is rougher on short length scales than one would 
expect from Eq. While the second observation is not 
unusual and has been predicted theoretically for sys- 
tems with short range and long range interactions j25j . 
the first is clearly characteristic for liquid crystal inter- 
faces. In Sec. 13 we will analyze it within the Landau-de 
Gennes theory. In particular, we will discuss the influence 
of elastic interactions. We find that the anisotropy of the 
spectrum can already been explained within an approxi- 
mation that excludes elastic interactions. However, adding 
the latter changes the spectrum qualitatively to the effect 
that the leading surface tension term becomes isotropic, 
and the anisotropy is governed by additional higher order 
terms. We summarize and conclude in Sec. 0] 



2 Berreman anchoring on rough and 
patterned surfaces 

We consider a nematic liquid crystal confined by a surface 
at z = h{x,y), which locally favors planar anchoring (i.e., 
alignment parallel to the surface) . The surface fluctuations 
h{x, y) are assumed to be small. The bulk free energy is 



given by the Frank elastic energy |Tll44j 

Ff = - J dxdy J dz i^Ki{Vnf + K2{n{V x n))^ 

+ if3(nx (Vxn))2}, (2) 

where n is the director, a vector of length unity which 
describes the local orientation of the liquid crystal, and Ki 
are the elastic constants (splay, twist and bend) . Since the 
surface favors planar alignment and the bulk fluctuations 
are small, we assume that the orientation of the director 
deep in the bulk, nf,, lies in the (x, y)-plane and that local 
deviations from n;, are small. Without loss of generality, 
we take nf, to point in the y direction. Hence we rewrite 
the director as 

n = {u, \/ \ — V? — ,v) ^ (3) 

and expand the free energy |(2J) up to second order in pow- 
ers of It, w, and h. This gives 

Fp J dxdy J dz ^Ki{d^u + d^vf 

+ K2{d,u - d^vf + Ksiidyuf + (dyvf)}. (4) 

Next we perform a two dimensional Fourier transform 
(x, y) — > q. Minimizing Fp in the bulk leads to the Euler- 
Lagrange equations 

( K2d,,-Kiql-K^ql -iq^{K2-Ki)d, \fu\ 
\ -iqAK2-Ki)d,_ K^d,,~K2ql-K3qlJ {v J 

For the boundary conditions {u,v) — > for z — > — oo and 
{u, v) = {uo, vq) at z = 0, the solution has the form 

('^\ ^ (^1^ -A2\ f ciexp(Aiz)\ 
\v J \\i iqx J \c2exp{\2z) J 

with the coefficients 

(""A^ I ( il^ >^2\ (uo\ .g. 

\C2 J X1X2 - q"^ ^Ix J \Vo J 

and the inverse decay lengths 

>^l2-ll + ll^- (7) 
Inserting that into the Frank energy one obtains 

^^ = 1 Aq /"^' 2 |Ai|^oP+ X2\vo\^+ 2q.x^{v*oUo)}, 

(8) 

where uo(q) and vo{q) are the values of u{q),v{q) at the 
surface. This is a general result, which we shall also use in 
Sec.El 

Now we study more specifically a liquid crystal in con- 
tact with a fixed patterned surface (fixed h{x,y)), which 
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anchors in an unspecifically planar way. The surface en- 
ergy is taken to be of Rapini Papoular type 01] 

F.^ao JdA{l + ^{noN)^) with > 0, (9) 



where dA = dx dy a/1 + (dxh)-^ + [dyhY is the local sur- 
face area element at {x,y), and 



N = 



(-dxK^dyh,!) (10) 



the local surface normal. Planar anchoring implies ag > 
0. As before, we rewrite the director at the surface no 
in terms of local deviations uq, vq, according to Eq. 
perform a Fourier transform {x^y) — > q, and expand Fs 
up to second order in uq, "Wo, and h. Omitting the constant 
contribution ctq^, this gives 



Fs = 



y"dq{|/i|\2^ao|i'o- jga^l'}. (11) 



We combine ^ and pi|) and minimize the total free en- 
ergy F = Fp + Fs with respect to uq and vq- The result 
is 



F^^Jdq\hfq^ {ao + Ksq- 



-2^ / -2\ 



with (iy — qy/q and 



k{q^^ = ll^l + qliKz/K.-l). 



(12) 



(13) 



The result can be generalized easily for the case that the 
bulk director points in an arbitrary planar direction 



rifc = (cos (/)o, sin 00,0) 



(14) 



by simply replacing qy with n;,q/g. 

Eq. (|12|l already shows that the bulk director will favor 
orientations where the amplitudes |/i(q)| arc small, i.e., 
the roughness is low. To quantify this further, we expand 
the integrand of (|12|l for small wave vectors in powers of q. 
The angle dependent part of the free energy as a function 
of the bulk director angle (pQ then takes the form 

^(0o) = ^/ d<j>cos^{(j)-(j)o)k{cos''{(j)-(j>o))H{<j>), (15) 



where the roughness spectrum |ft,(q) ^ enters the anchoring 
energy solely through the function 



H{cf)^ / dqq^\h{ci)\' 
Jo 



(16) 



It is convenient to expand H{(f>) into a Fourier series with 
coefficients 



Hn — 



1 

2^ 



,iy(0)e-'"'^ -|i/„|e™"''. (17) 



Similarly, we define 

Cn = — [dip cosV ^;(cos^(/)) 
2vr J 



(18) 



The coefficients c„ are real and vanish for odd n. In the 
case K3 = Ki {e.g., in the Landau-de Gennes theory. Sec. 
13.1(1 . one has k = 1, and the series c„ stops at \n\ =4 with 
C2 = 1/4 and C4 = 1/16. In real materials PP, the elastic 
constant K3 is typically larger than Ki by a factor of 1-3, 
and the series does not stop. However, the coefficients for 
\n\ < 4 remain positive, and the coefficients for \n\ > 4 
become very small, such that they may be neglected. 

Omitting constant terms that do not depend on (j)Q, 
the anchoring energy can then be written as 



F{M = -TT-fCs ^ Cn\H„ \cosn{<j)o - a„). 



(19) 



The anchoring angle is the angle that minimizes F{(f)Q). 
In general, the n = 2 term will dominate, and one gets 
approximately 



a2 = -ar^ 



(20) 



We note that the angles a„ in Eq. 1)17(1 correspond to di- 
rections where the height fluctuations are small, because 
the contributions of Hq and Hn to the spectral function 
-ff(0) have opposite signs. Hence Eq. (pn|l impHes that 
the surface aligns the nematic fluid in a direction where 
the surface is smooth. At given anchoring angle (po, we 
can also calculate the anchoring strength. To this end, 
we expand the anchoring energy about (/>o and obtain 
F(0o) — F(0o) + -^(0o — ^0)^ with the anchoring strength 



n=2.4 



\H„\cosn{ 



0- (21) 



We conclude that elastic interactions in nematic liquid 
crystals on anisotropically rough surfaces induce an an- 
choring energy in a direction of low roughness. The cen- 
tral quantities characterizing the surface roughness are the 
two coefficients _ff2,4 defined by Eqs. ((16(1 and l(17() . These 
quantities determine the anchoring strength (Eq. 1(2111 ). 
and the anchoring angle (Eq. 1(20(1 ). The anchoring mech- 
anism only requires an unspecific tendency of the liquid 
crystal to align parallel to the interface (Eq. lO). Given 
such a tendency, the anchoring energy no longer depends 
on the surface parameters, ao and uq. The only relevant 
material parameters are the splay and bend elastic con- 
stants in the bulk, Ki and K^, and the squared surface 
anisotropy, which is characterized by the coefficients iJ2,4. 

We note that our treatment premises that the nematic 
liquid stays perfectly ordered at the surface. In reality, 
rough surfaces may reduce the order, which in turn influ- 
ences the anchoring properties |46II47[ . This has not been 
considered here. 
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3 Capillary waves at the nematic/isotropic 
interface 

In this section we study the capillary wave spectrum of 
freely undulating nematic/isotropic (NI) interfaces. The 
problem is similar to that considered in the previous sec- 
tion (Sec. HJ, with two differences: (i) The interface po- 
sition h{x, y) is free and subject to thermal fluctuations, 
and (ii) the nematic order at the interface drops smoothly 
to zero. The second point implies, among other, that the 
elastic constants are reduced in the vicinity of the inter- 
face. 

In many systems, the anchoring at Nl-interfaces is pla- 
nar. As a zeroth order approach, we neglect the softness of 
the profile and approximate the interfacial structure by a 
steplike structure (sharp-kink approximation), and the in- 
terfacial free energy by Eq. (|12|l with effective parameters 
(To and ao- Generally the capillary waves of an interface 
with an interfacial free energy of the form 



F = \j dqlMq)P^(q) 

are distributed according to 

(|Mq)P>-fcBT/i:(q)). 



(22) 



(23) 

Thus the free energy (|12|l yields the capillary wave spec- 
trum 



(IMq)P) 



2 , 3 -^3 .4 - 

q +q qyK 

CTO ^ 



(24) 



As before, qy is the component of the unit vector q/q in 
the direction of the bulk director. 

The result (|24|l shows already the three remarkable fea- 
tures, which will turn out to be characteristic for the NI 
interface. First, the capillary wave spectrum is anisotropic, 
the capillary waves in the direction parallel to the direc- 
tor (^j^) being smaller than in the direction perpendicular 
to the director. Second, the leading (quadratic) term is 
still isotropic; the anisotropy enters through the higher 
order terms. Third, in contrast to simple fluids with short 
range interactions, the capillary wave spectrum cannot be 
expanded in even powers of but it contains additional 
cubic (and higher order odd) terms. This implies that the 
capillary wave spectrum is nonanalytic in the limit q ^ 0. 

These findings are gratifying. However, the sharp kink 
description of the NI interface is inadequate. The Frank 
free energy ^ describes nematic liquid crystals with con- 
stant local order parameter, whereas at NI interfaces, the 
nematic order parameter drops softly to zero. Moreover, 
the surface anchoring at NI interfaces is an intrinsic prop- 
erty of the interface, which depends itself on the local elas- 
tic constants. We will now consider our problem within a 
unified theory for nematic and isotropic liquid crystals, 
the Landau-de Gennes theory. 

3.1 Landau-de Gennes theory 

The Landau-de Gennes theory is based on a free energy 
expansion in powers of a symmetric and traceless (3 x 3) 



order tensor field Q(r). 

F = |dr{^Tr(Q^) + |Tr(Q^)-H^Tr(Q^)^4-^Tr(Q^) 

+ ^drQjkdiQjk + ^diQijdkQkj^ (25) 

Following a common assumption, we neglect the possibil- 
ity of biaxiality and rewrite the order tensor as 0S1 



?^j(i") = \s{Y){2,n,{Y)nj{Y) - (5y). 



(26) 



Here S is the local order parameter, and n a unit vec- 
tor characterizing the local director. In the homogeneous 
case {diQjk = 0), the free energy H25|l predicts a first or- 
der transition between an isotropic phase (I) with S = 
and an oriented nematic phase (N) with S — Sq — 
-2/9B/(Ci + C2/2). We recall briefly the properties of 
a flat NI interface at coexistence for a system with fixed 
director n, as obtained from minimizing H25(l |48| : The 
order parameter profile has a simple tanh form 

S{z) = SqS{z/0 with 5(r) ^ l/(e^ + 1). (27) 

The interfacial width 



and the interfacial tension 



a = crov/l + "(nN)2 



3(Ci + C2/2) 4 , , - 
0^0 ^0^0 I (29) 



both depend in the same way on the angle between the 
director n and the surface normal N, via the parameter 



1 



2 Li 



L2/6 



(30) 



The quantity (Tq sets the energy scale, ^0 the length 
scale, and the "order parameter scale". (Note that S 
can be rescaled even though it is dimensionless) . Hence 
only one characteristic material parameter remains, e.g., 
the parameter a. In the following, we shall always use 
rescaled quantities S ^ S/Sq, length — » length/^O; energy - 
energy/tTo- The free energy at coexistence can then be 
rewritten as EHI 



j dv{f 



+ 91 +92+ 93 + .94} 



(31) 



with / = 3S'^(S'^ - 1) (at coexistence) 

51 = 3((V5)2+a(nV5)2) 

52 = 12a((Vn)(nVS') i(n x V x n)(VS)) 

53 = 3(^(3 + 2a)(Vn)2 + (3 - a)(n ■ V x n)^ 

-t- (3 + 2a)(n x V x n)^^ 
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The first term f{S) describes the bulk coexistence, the 
middle terms gi and 172 determine the structure of the 
interface, and the last term establishes the relation to the 
Frank elastic energy, Eq. ^ . We note that in this version 
of the Landau-de Gennes theory, the splay and the bend 
elastic constants are identical, Ki — K3, hence k{qy) = 1 
in Eq. if^ . 

Eq. H3f|l will be our starting point. As in Sec. El we 
will assume without loss of generality that the interface is 
on average located at z = 0, and that the bulk director 
far from the surface points in the y-direction. 



3.2 Constant director approximation 



As before, we assume that the profile has the form S = 
S((z — h)/^). After inserting this Ansatz in (|31() . expand- 
ing in 9 up to second order, and minimizing with respect 
to 5, a lengthy calculation yields the surface free energy 
(omitting constant terms) 



F 



\ [dci\h\'{c 



3c' 



ci + C2q: 



y 



esq; 



:)}(35) 



with 



1 + K 



- 2k 



Cl 



C2 



3 

f + 2k 
3 + 2a / 1 



2k 



, 3 -I- a / 1 



K 



3a V2k-1 



r 
C3 



a \K — f 

3-a / 1 

" ,2k - 1 



a 



We return to considering undulating interfaces with vary- 
ing position h{x, y). In the simplest approach, the director 
is constant throughout the system, n = (0,1,0). Elastic 
interactions are thus disregar_ded. For the order parameter, where we have defined generalized binomial coefficients 
we make the Ansatz S{v) = S{{z — h{x,y)) / , where S is 
the tanh profile from Eq. (|27|l . and the interfacial width 
^ varies with the local surface normal N H10() according 
to Eq. H28|l . Inserting this into the free energy (|31|l . and 
Fourier transforming {x,y) — > q, one obtains 



r{n+l) 

r{a + i)r{n- a + iy 



(32) 



This result is quite robust. As we shall see in Sec. H3.4|l . 
it is also obtained with a rather different Ansatz for S', as 
long as the director is kept constant. We can now apply 
Eq. (ESI and obt ain the capillary wave spectrum 



(IMq)l') 



(g2 + aql) 



(33) 



As a consistency check, we also inspected the result for 
9{x, y, z) directly. It is proportional to dyh as expected. 

The comparison of the free energy H35(l with the corre- 
sponding result for fixed director, Eq. (|32|) . shows that the 
anisotropy of the surface fluctuations is reduced. For a fur- 
ther analysis, it would be necessary to minimize the free 
energy expression with respect to the variational pa- 
rameter K. Unfortunately, k enters in such a complicated 
way, that this turns out to be unfeasible. Numerically, we 
find that the capillary wave spectrum, obtained via Eq. 
(I23|l . varies only little with k. For any reasonable value 



of K, 



> 2, the result differs from that obtained 



(recalling that wave vectors q are given in units of l/^o 
H28(l '). which is anisotropic. 

Hence already this simple approach predicts anisotropic 
capillary wave amplitudes. The capillary waves are weak- 
est in the y direction, which is the direction of the bulk 
director. It is interesting that we get the anisotropy al- 
ready at this point. One might have suspected that the 
dampening of waves parallel to the director is caused by 
director-director interactions. This turns out not to be 
the case, instead an interaction between the director and 
the order parameter gradient nSI S is responsible for the 
anisotropy. As the director wants to align parallel to the 
surface, waves parallel to the director have higher energy. 



3.3 Relaxing the director: A variational approach 

Since the interface locally favors parallel anchoring, one 
would expect that the director follows the undulations of 
the interface (Fig. ^ . This motivates a variational Ansatz 
n — (O,cos0,sin0) with 



e{x,y,z) 



g{x,y) exp[-(z 



h{x,y))]. (34) 



with the constant director approximation (Eq. H33() l by 
less than one percent. Within the present approximation, 
the effect of relaxing the director is negligeable. This is 
mostly due to the fact that Eq. H34(l still imposes rather 
rigid constraints on the director variations in the nematic 
fluid. 



3.4 Local profile approximation 

A more general solution can be obtained with the addi- 
tional approximation that the width of the interface is 
small, compared to the relevant length scales of the inter- 
facial undulations. In that case, the interface and the bulk 
can be considered separately, and we can derive analytical 
expressions for the capillary wave spectrum. The assump- 
tion that length scales can be separated is highly ques- 
tionable, because interfacial undulations are present on 
all length scales down to the molecular size. Nevertheless, 
computer simulations of other systems (Ising models 00] 
and polymer blends |36| ) have shown that intrinsic profile 
models can often describe the structure of interfaces quite 
successfully. 

We separate the free energy (|31|) into an interface and 
bulk contribution, F = Fs + Fp. The bulk contribution 
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/ N 




gis 


r=(x,y,z) 








Nematic 


h(x',y') - 


— > 


923 



Fig. 1. Nematic-Isotropic interface with local coordinates. 



Fp has the form (|2Jl with the elastic constants Ki = = 
6(3 + 2a) and K2 = 6(3 — a), and accounts for the elastic 
energy within the nematic region. The integrand in the 
expression for the remaining surface free energy Fs van- 
ishes far from the surface. We assume that the focal order 
parameter profile has mean-field shape in the direction 
perpendicular to the surface. More precisely, we make the 
Ansatz 



s{v) = 5(c/C), 



vs* 



1 dS 



N, (36) 



(cf. H27|) . H28|l 'l. where N is the focal surface normal as 
usual, and ( is the distance between r and the closest inter- 
face point. The {x, y) coordinates at this point are denoted 
(a;', y'). (see Fig.P). To evaluate Fs, we make a coordinate 
transformation r — > {x' ,y', and integrate over The re- 
lation between the coordinates is r = {x' ,y' , h{x' ,y')) + 
NC, and the Jacobi determinant for the integral is in sec- 
ond order of h 

^+]^{{dxhf +{dyhY)-C{dxxh+dyyh)+C,'^ [hxxhyy-hxyhyx) 

We begin with reconsidering the constant director case, 
n = const. The Frank free energy then vanishes, and the 
surface free energy takes exactly the form of Eq. 
Hence the present approximation leads to the same ex- 
pression as the approximation taken in Sec. 13.21 This un- 
derlines the robustness of the result H33|l . 

In the general case, we must make an Ansatz for the 
variation of the director n in the vicinity of the surface. 
We assume that it varies sufficiently slowly, so that we can 
make a linear approximation 



n(x',y',C)«n(a;',y',0)-f C(N-V)n. 



(37) 



As before, we take the bulk director to point in the y- 
direction. The local director deviations in the x and z 
direction are parametrized in terms of two parameters u 
and V according to Eq. @. 

After inserting Eqs. and ||27Il mto Eq. (PIJ, ex- 
panding up to second order in h, u, and v, and some par- 
tial integrations, we obtain the surface free energy 



(38) 



Fs ^ dx dy {gis + 52s + 53s) 



7r2 -6 



-{dxh){dyUQ) - void^uo - Sd^uo) 
+2{d,vo){dxUQ + d^vo) + 2{dzUQ){dzUQ - dxVo)^ 

gas = s(i3 + 2a){dxUo + d,VQf + {3-a){dzUo-dxVo)'^ 

+ {3 + 2a){{dyUof + {dyVo'f)] 



where Uq and ug are the values of u and v at the inter- 
face. The first contribution gis describes the effect of the 
anisotropic local surface tension. The second contribution 
(72s arises from the coupling between the director varia- 
tions diU with the order parameter variation \/S at the 
interface. The last contribution g^s accounts for the re- 
duction of the Frank free energy in the interface region. 

The ensuing procedure is similar to that of Sec. El We 
first minimize the bulk free energy Fp, which leads in the 
most general case to Eq. ()12(l . Then we minimize the total 
energy F = Fp + Fs with respect to uq and wq, using 
Eqs. (O and ^ to estimate the derivatives dzUo, dzVQ. 
The result has the form (|22() and gives the capillary wave 
spectrum via (|23l) . Unfortunately, the final expression is a 
rather lengthy, and we cannot give the formula here. We 
will discuss it further below. 

A more concise and qualitatively similar result is ob- 
tained with the additional approximation, that the direc- 
tor only varies in the z direction, u = 0. In that case, the 
minimization of the Frank free energy with respect to 
V yields v{q, z) = vq cxp(gA:z) with 



V 3 + 2a 



1 



-2 3a 
Qy: 



'3-a 



(39) 



(My ~ at 2: = we have 9zUo — qkuQ. The 

Frank energy takes the form 

1 



Fp^ — I dq{3 + 2a)qkiq'y)\vo\ 
12a 



(40) 



This equation replaces Eq. H12|) . The surface free energy 
Fs H38|l is also greatly simplified. After minimizing the 
sum F = Fp + Fs with respect to vq and applying Eq. 
(|23|l . one obtains the capillary wave spectrum 



ksT /(To 2 , 2 
- q +aqy 



(IMq)P) 

xfl- 



(1 - Qqk - 'iq^f 



l + 185fc((i +i)(l_2gfc) + ^gfc) 
Expanding this for small wavevectors q gives 
ksT/cTQ 



(41) 



(IMq)l^) 



g' + C/g' 18(l + a)A:-g*---), (42) 
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(q g 

Fig. 2. Capillary wave spectrum in the direction perpendic- 
ular to the bulk director (qx) (solid line) and parallel to the 
bulk director (qy) as obtained from the constant director ap- 
proximation (dotted line) and the local profile approximation 
(dashed line). In the perpendicular direction, both approaches 
give the same result. Note the presence of an unphysical pole 
at (q^o)^ = 0.24. The inset shows a blowup for very small 
g-vectors, illustrating that the spectrum becomes isotropic in 
the limit g — > 0. The parameter is a = 3/13, corresponding to 
Li = 2L2. 



where the coefficient of the fourth order term is negative. 

Comparing this solution to Eq. , one notes obvious 
qualitative differences. In the constant director approxi- 
mation, Eq. H33|l . one has an anisotropic effective surface 
tension: The capillary wave spectrum has the form Q 
with a — <7o{l + aQy). The present treatment shows that 
the elastic interactions remove the anisotropy in the sur- 
face tension term (order q^), but introduce new anisotropic 
terms that are of higher order in q. This is consistent with 
the preliminary results from our earlier zeroth order ap- 
proach, Eq. f^ . 

We turn to discussing the full solution of the local pro- 
file approximation, where both variations of u and v are 
allowed. In the directions parallel and perpendicular to 
the bulk director (x and y), the capillary wave spectrum 
turns out to be the same as in (|^ . It is shown in Fig. El 
for a typical value of a (taken from Ref. |3H1) and com- 
pared to the constant director approximation. The capil- 
lary waves in the direction perpendicular to the bulk direc- 
tor (the x-direction) are not influenced by the elastic in- 
teractions: The amplitudes only contain a q^ contribution 
and are identical in the constant director approximation 
and the local profile approximation. In contrast, the cap- 
illary waves in the direction parallel to the bulk director 
(the y-direction) are to a large extent dominated by the 
cubic term. The fact that the spectrum becomes isotropic 
in the limit g — > becomes only apparent at very small q 
vectors, q£_o < 0.005 (see inset of Fig. 12)). 

The effect of relaxing both u and v in Eq. (PJ becomes 
apparent when examining directions that are intermediate 
between x and y. The expansion of the full solution in 
powers of q gives an expression that is similar to Eq. (|42|) , 



(IMq)P) 



= q'+ qliq' CM) + / C^iql) + • • • ) (43) 



but has different coefficients Fig. O shows the coeffi- 
cients C3 and C4 as a function of qy and compares them 
with the corresponding quantities obtained from the sim- 
plified solution 1)41(1 . In the full solution, the Ci stay much 
smaller in the vicinity oi qy 0. This demonstrates once 
more that elastic interactions reduce the anisotropy of the 
capillary waves. 
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-20 
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Fig. 3. Coefficients of (a) the cubic term and (b) the fourth 
order term vs. = qy/q in the expansion 1431 and 14211 . at 
a = 3/13. 



Unfortunately, we also find that the free energy F — 
Ff + Fs (from Eqs. (O and (|2H1)), as a functional of 
wo(q), i'oIq); becomes unstable for larger wavevectors q. 
The instability gives rise to the unphysical pole at (g^^o)^ — 
0.24 in Fig. 121 In the direction qy = 1, the pole is encoun- 
tered at g^o ~ 0-5 for all a. In other directions, it moves 
to even smaller q values, which reduces the stability fur- 
ther. At a = 3/13, the free energy is stable only up to 
9^0 ~ 0.023, corresponding to a length scale of ^ 270^o- 
The approximation is bound to break down on smaller 
length scales. 

Hence the region of validity of the theory is very small. 
The instability is presumably caused by the way the lo- 
cal profile approximation was implemented. In particular, 
our estimate for the values of the partial derivatives dzUo 
and dzVo at the interface must be questioned. They were 
obtained from extrapolating the bulk solution, which is 
however only valid for constant, saturated, order param- 
eter 5 = 5*0. In order to assess the effect of this con- 
straint, we have considered a second approximation: The 
derivative 9^11 at the interface is taken to be independent 
of the bulk solution. In the interface region, the director 
is assumed to vary linearly with the slope cJ^n. At the 
distance ^ 2.5^o from the interface, the proffie changes 
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continuously to the exponential bulk solution. In this ap- 
proximation, the value of the derivative c^^n at the surface 
is an additional variable, which can be optimized inde- 
pendently. We have minimized the free energy with this 
Ansatz and u = (i.e., no director variations in the x- 
direction), and calculated the capillary wave spectrum. 
The result is shown in Fig. 0] The instability from Fig. El 
disappears. The other characteristic features of the spec- 
trum remain. It becomes isotropic for very small wavevec- 
tors, q < 0.003^0, corresponding to length scales of several 
thousand correlation lengths ^o- On smaller length scales, 
it is anisotropic. It is worth noting that at least on this 
level of approximation, the capillary waves in the direc- 
tion perpendicular to the director (the x direction) are 
not affected by the elastic interactions. 



6 




12 3 4 

Fig. 4. Capillary wave spectrum from the local profile approx- 
imation with independent surface derivative dzU, compared to 
the constant director approximation, in the directions parallel 
and perpendicular to the bulk director. The pole in Fig.|5|has 
disappeared. At small wavevectors, the curves are the same as 
in the approximation of Fig. |5| 



4 Summary and discussion 

To summarize, we have studied the interplay of elastic 
interactions and surface undulations for nematic liquid 
crystals at rough and fluctuating interfaces using appro- 
priate continuum theories: the Frank elastic energy and 
the Landau-de Gennes theory. 

In the first part, we have considered nematic liquid 
crystals in contact with a surface of given geometry, char- 
acterized by a fixed height function h{x,y). We have re- 
analyzed the effect of Berreman anchoring, i.e., the phe- 
nomenon that elastic interactions are sufficient to align 
liquid crystals exposed to anisotropically rough surfaces. 
Our treatment allowed to derive explicit equations for 
the anchoring angle and the anchoring strength for given 
(arbitrary) height function h{x,y). In particular, we find 
that the resulting azimuthal anchoring coefficient depends 



only on the surface anisotropy and the bulk elastic con- 
stants Ki and K^, and not on the chemical surface inter- 
action parameters such as the interfacial tension and the 
zenithal anchoring strength. The contribution of the sur- 
face anisotropy to the anchoring energy has recently been 
verified by Kumar et al • We hope that our results will 
stimulate systematic experimental research on the role of 
the elastic constants as well. 

In the second part, we have examined the inverse prob- 
lem, the effect of the nematic order on capillary wave fiuc- 
tuations of NI interfaces. The work was motivated by a 
previous simulation study, where it was found that the 
capillary wave amplitudes were different in the direction 
parallel and perpendicular to the bulk director. Our anal- 
ysis shows that this effect can be understood within the 
Landau-de Gennes theory. As in the simulation, the waves 
parallel to the director are smaller than those perpendicu- 
lar. The anisotropy is caused by a coupling term between 
the director and the order parameter gradient (nVS), 
which locally encourages the director to align parallel to 
the surface, and thus penalizes interfacial undulations in 
the direction of the director. 

The infiuence of elastic interactions mediated by the 
nematic bulk fluid was investigated separately. We find 
that they reduce the anisotropy and change the capillary 
wave spectrum qualitatively. In the absence of elastic in- 
teractions, i.e., with fixed director, the anisotropy man- 
ifests itself in an anisotropic surface tension. If one al- 
lows for director modulations, the surface tension becomes 
isotropic, and the anisotropy is incorporated in higher or- 
der terms in the wave vector q. In particular, we obtain 
a large anisotropic cubic term. The fourth order term is 
generally negative, i.e., we have a negative "bending rigid- 
ity" . This is consistent with previous observations from 
simulations As we have noted in the introduction, 
bending rigidities are often negative in fiuid-fluid inter- 
faces, for various reasons. In the case of nematic/isotropic 
interfaces, the elasticity of the adjacent medium provides 
an additional reason. 

We have shown that the higher order terms dominate 
the capillary wave spectrum on length scales up to several 
thousand correlation lengths. This has serious practical 
consequences. The analysis of capillary waves is usually a 
valuable tool to determine interfacial tensions from com- 
puter simulations. In liquid crystals, however, this method 
must be applied with caution. More generally, our re- 
sult suggests that the apparent interfacial tension of ne- 
matic/isotropic interfaces should be strongly affected by 
finite size effects. In fact, Vink and Schilling have recently 
reported that the interfacial tension obtained from com- 
puter simulations of soft spherocylinders varies consider- 
ably with the system size '5D"^r . 

We thank Marcus Miiller for a useful comment, and 
the German Science Foundation for partial support. 
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